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Abstract Analogue models of gravity have provided an experimentally realizable 
0^ . test field for our ideas on quantum field theory in curved spacetimes but they have 

also inspired the investigation of possible departures from exact Lorentz invariance 
at microscopic scales. In this role they have joined, and sometime anticipated, sev- 
O ' eral quantum gravity models characterized by Lorentz breaking phenomenology. A 

O"^ crucial difference between these speculations and other ones associated to quan- 

!— i ' turn gravity scenarios, is the possibility to carry out observational and experimen- 

tal tests which have nowadays led to a broad range of constraints on departures 
from Lorentz invariance. We shall review here the effective field theory approach to 
C*~) ', Lorentz breaking in the matter sector, present the constraints provided by the avail- 

*> • able observations and finally discuss the implications of the persisting uncertainty 

on the composition of the ultra high energy cosmic rays for the constraints on the 
higher order, analogue gravity inspired, Lorentz violations. 



1 Introduction 



Our understanding of the fundamental laws of Nature is based at present on two 
different theories: the Standard Model of Fundamental Interactions (SM), and clas- 
rS ■ sical General Relativity (GR). However, in spite of their phenomenological success, 

jrt , SM and GR leave many theoretical questions still unanswered. First of all, since we 

feel that our understanding of the fundamental laws of Nature is deeper (and more 
accomplished) if we are able to reduce the number of degrees of freedom and cou- 
pling constants we need in order to describe it, many physicists have been trying to 
construct unified theories in which not only sub-nuclear forces are seen as different 
aspects of a unique interaction, but also gravity is included in a consistent manner. 
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Another important reason why we seek for a new theory of gravity comes di- 
rectly from the gravity side. We know that GR fails to be a predictive theory in some 
regimes. Indeed, some solutions of Einstein's equations are known to be singular at 
some points, meaning that in these points GR is not able to make any prediction. 
Moreover, there are apparently honest solutions of GR equations predicting the ex- 
istence of time-like closed curves, which would imply the possibility of traveling 
back and forth in time with the related causality paradoxes. Finally, the problem of 
black-hole evaporation seems to clash with Quantum Mechanical unitary evolution. 

This long list of puzzles spurred an intense research toward a quantum theory of 
gravity that started almost immediately after Einstein's proposal of GR and which 
is still going on nowadays. The quantum gravity problem is not only conceptually 
challenging, it has also been an almost metaphysical pursue for several decades. 
Indeed, we expect QG effects at experimentally/observationally accessible energies 
to be extremely small, due to suppression by the Planck scale M v \ = y/Tic/Gn ~ 
1.22 x 10 19 GeV/c 2 . In this sense it has been considered (and it is still considered 
by many) that only ultra-high-precision (or Planck scale energy) experiments would 
be able to test quantum gravity models. 

It was however realized (mainly over the course of the past decade) that the sit- 
uation is not as bleak as it appears. In fact, models of gravitation beyond GR and 
models of QG have shown that there can be several of what we term low energy 
"relic signatures" of these models, which would lead to deviation from the standard 
theory predictions (SM plus GR) in specific regimes. Some of these new phenom- 
ena, which comprise what is often termed "QG phenomenology", include: 

• Quantum decoherence and state collapse 0] 

• QG imprint on initial cosmological perturbations J2] 

• Cosmological variation of couplings |3]|U 

• TeV Black Holes, related to extra-dimensions J5] 

• Violation of discrete symmetries J6) 

• Violation of space-time symmetries Q 

In this lecture I will focus upon the phenomenology of violations of space-time 
symmetries, and in particular of Local Lorentz invariance (LLI), a pillar both of 
quantum field theory as well as GR (LLI is a crucial part of the Einstein Equivalence 
Principle on which metric theories of gravity are based). 



2 A brief history of an heresy 

Contrary to the common trust, ideas about the possible breakdown of LLI have 
a long standing history. It is however undeniable that the last twenty years have 
witnessed a striking acceleration in the development both of theoretical ideas as 
well as of phenomenological tests before unconceivable. We shall here present an 
incomplete review of these developments. 
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2.1 The dark ages 

The possibility that Lorentz invariance violation (LV) could play a role again in 
physics dates back by at least sixty years [8] |9] \W\ QT| [T2] Q~3) and in the seventies 
and eighties there was already a well established literature investigating the possible 
phenomenological consequences of LV (see e.g. |[T^[l31[T6l[T7l[T8l[T9l ). 

The relative scarcity of these studies in the field was due to the general expecta- 
tion that new effects were only to appear in particle interactions at energies of order 
the Planck mass M p \. However, it was only in the nineties that it was clearly realized 
that there are special situations in which new effects could manifest also at lower 
energy. These situations were termed "Windows on Quantum Gravity". 



2.2 Windows on Quantum Gravity 

At first glance, it appears hopeless to search for effects suppressed by the Planck 
scale. Even the most energetic particles ever detected (Ultra High Energy Cosmic 
Rays, see, e.g., fl20ll2Tl )haveE' < 10 11 GeV ~ 10~ 8 M p /. However, even tiny correc- 
tions can be magnified into a significant effect when dealing with high energies (but 
still well below the Planck scale), long distances of signal propagation, or peculiar 
reactions (see, e.g., |(7) for an extensive review). 
A partial list of these windows on QG includes: 

• sidereal variation of LV couplings as the lab moves with respect to a preferred 
frame or direction 

• cumulative effects: long baseline dispersion and vacuum birefringence (e.g. of 
signals from gamma ray bursts, active galactic nuclei, pulsars) 

• anomalous (normally forbidden) threshold reactions allowed by LV terms (e.g. pho- 
ton decay, vacuum Cherenkov effect) 

• shifting of existing threshold reactions (e.g. photon annihilation from Blazars, 
ultra high energy protons pion production) 

• LV induced decays not characterized by a threshold (e.g. decay of a particle from 
one helicity to the other or photon splitting) 

• maximum velocity (e.g. synchrotron peak from supernova remnants) 

• dynamical effects of LV background fields (e.g. gravitational coupling and addi- 
tional wave modes) 

It is difficult to assign a definitive "paternity" to a field, and the so called Quan- 
tum Gravity Phenomenology is no exception in this sense. However, among the pa- 
pers commonly accepted as seminal we can cite the one by Kostelecky and Samuel 
ll22ll that already in 1989 envisaged, within a string field theory framework, the pos- 
sibility of non-zero vacuum expectation values (VEV) for some Lorentz breaking 
operators. This work led later on to the development of a systematic extension of the 
SM (what was later on called "minimal standard model extension" (mSME)) incor- 
porating all possible Lorentz breaking, power counting renormalizable, operators 



Stefano Liberati 



(i.e. of mass dimension < 4), as proposed by Colladay and Kostelecky 1231 . This 
provided a framework for computing in effective field theory the observable conse- 
quences for many experiments and led to much experimental work setting limits on 
the LV parameters in the Lagrangian (see e.g. 11241 for a review). 

Another seminal paper was that of Amelino-Camelia and collaborators ll25l 
which highlighted the possibility to cast observational constraints on high energy 
violations of Lorentz invariance in the photon dispersion relation by using the afore- 
mentioned propagation over cosmological distance of light from remote astrophys- 
ical sources like gamma ray bursters (GRBs) and active galactic nuclei (AGN). The 
field of phenomenological constraints on quantum gravity induced LV was born. 

Finally, we should also mention the influential papers by Coleman and Glashow 
l26ll27l[28l which brought the subject of systematic tests of Lorentz violation to the 
attention of of the broader community of particle physicists. 

Let me stress that this is necessarily an incomplete account of the literature which 
somehow pointed a spotlight on the investigation of departures from Special Relativ- 
ity. Several papers appeared in the same period and some of them anticipated many 
important results, see e.g. ll29l l30l . which however at the time of their appearance 
were hardly noticed (and seen by many as too "exotic"). 

In the years 2000 the field reached a concrete maturity and many papers pursued 
a systematization both of the framework as well as of the available constraints (see 
e.g. BT1 l32l [33ll ). In this sense another crucial contribution was the development 
of an effective field theory approach also for higher order (mass dimension greater 
than four), naively non-power counting renormalizable, operatorsLj. This was firstly 
done for dimension 5 operators in QED 11371 by Myers and Pospelov and later on 
extended to dimension 6 operators by Mattingly ll38l . 

Why all this attention to Lorentz breaking tests developed in the late nineties and 
in the first decade of the new century? I think that the answer is twofold as it is re- 
lated to important developments coming from experiments and observation as well 
as from theoretical investigations. It is a fact that the zoo of quantum gravity mod- 
els/scenarios with a low energy phenomenology had a rapid growth during those 
years. This happened mainly under the powerful push of novel puzzling observa- 
tions that seemed to call for new physics possibly of gravitational origin. For ex- 
ample, in cosmology these are the years of the striking realization that our universe 
is undergoing an accelerated expansion phase [39. 40| which apparently requires a 
new exotic cosmological fluid, called dark energy, which violates the strong energy 
condition (to be added to the already well known, and still mysterious, dark matter 
component). 

Also in the same period high energy astrophysics provided some new puzzles, 
first with the apparent absence of the Greisen-Zatsepin-Kuzmin (GZK) cut off f4~Tl 
l42l (a suppression of the high-energy tail of the UHECR spectrum due to UHECR 
interaction with CMB photons) as claimed by the Japanese experiment AGASA 



1 Anisotropic scaling |34||35|[36l techniques were recently recognized to be the most appropriate 
way of handling higher order operators in Lorentz breaking theories and in this case the highest 
order operators are indeed crucial in making the theory power counting renormalizable. This is 
why we shall adopt sometime the expression "naively non renormalizable" 
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114311 , later on via the so called TeV-gamma rays crisis, i.e. the apparent detection of 
a reduced absorption of TeV gamma rays emitted by AGN [44 j. Both these "crises" 
later on subsided or at least alternative, more orthodox, explanations for them were 
advanced. However, their undoubtedly boosted the research in the field at that time. 

It is perhaps this past "training" that made several exponents of the quantum 
gravity phenomenology community the among most ready to stress the apparent 
incompatibility of the recent CERN-LNGS based experiment OPERA B31 measure 
of superluminal propagation of muonic neutrinos and Lorentz EFT (see e.g. l46ll47l 
|48]|49)- There is now evidence that the Opera measurement might be flawed due to 
unaccounted experimental errors and furthermore it seems to be refuted by a similar 
measurement of the ICARUS collaboration ll50l . Nonetheless, this claim propelled a 
new burst of activity in Lorentz breaking phenomenology which might still provide 
useful insights for future searches. 

Parallel to these exciting developments on the experimental/observational side, 
also theoretical investigations provided new motivations for Lorentz breaking searches 
and constraints. Indeed, specific hints of LV arose from various approaches to Quan- 
tum Gravity. Among the many examples are the above mentioned string theory 
tensor VEVs l22l and space-time foam models llBTl l25l l52l l53l l54l . then semi- 
classical spin-network calculations in Loop QG 11551 . non-commutative geome- 
try J56 57 58 1, some brane-world backgrounds [59|. 

Indeed, during the last decades there were several attempts to formulate alterna- 
tive theories of gravitation incorporating some form of Lorentz breaking, from early 
studies ll60l I6T1 l62l l63l l64l to full-fledged theories such as the Einstein- Aether the- 
ory l65l l66l l67l and Hof ava-Lifshitz gravity ll3~5l l68l 1691 (which in some limit can 
be seen as a UV completion of the Einstein- Aether framework lf70lD . 

Finally, a relevant part of this story is related to the vigorous development in the 
same years of the so called condensed matter analogues of "emergent gravity" irTTl 
which is the main topic of this school. Let us then consider these models in some 
detail and discuss some lesson that can be drawn from them. 



3 Bose-Einstein condensates as an example of emergent Local 
Lorentz invariance 

Analogue models for gravity have provided a powerful tool for testing (at least in 
principle) kinematical features of classical and quantum field theories in curved 
spacetimes l7T1l . The typical setting is the one of sound waves propagating in a per- 
fect fluid Il72ll73l . Under certain conditions, their equation can be put in the form of 
a Klein-Gordon equation for a massless particle in curved spacetime, whose geome- 
try is specified by the acoustic metric. Among the various condensed matter systems 
so far considered, Bose-Einstein condensate (BEC) ll74l l75l had in recent years a 
prominent role for their simplicity as well as for the high degree of sophistication 
achieved by current experiments. In a BEC system one can consider explicitly the 
quantum field theory of the quasi-particles (or phonons), the massless excitations 
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over the condensate state, propagating over the condensate as the analogue of a 
quantum field theory of a scalar field propagating over a curved effective space- 
time described by the acoustic metric. It provides therefore a natural framework 
to explore different aspects of quantum field theory in various interesting curved 
backgrounds (for example quantum aspects of black hole physics l76l 1771 or the 
analogue of the creation of cosmological perturbations ll78l 1791 [80l ISTl ) or even, 
and more relevantly for our discussion here, emerging space time scenarios. 

In BEC, the effective emerging metric depends on the properties of the conden- 
sate wave-function. One can expect therefore the gravitational degrees of freedom 
to be encoded in the variables describing the condensate wave-function [751 , which 
is solution of the well known Bogoliubov-de Gennes (BdG) equation. The dynam- 
ics of gravitational degrees of freedom should then be inferred from this equation, 
which is essentially non-relativistic. The "emerging matter", the quasi-particles, in 
the standard BEC, are phonons, i.e. massless excitations described at low energies 
by a relativistic (we shall see in which sense) wave equation, however, at high en- 
ergies, the emergent nature of the underlying spacetime becomes evident and the 
relativistic structure of the equation broken. Let's see this in more detail as a con- 
ceptual exercise and for highlighting the inspirational role played in this sense by 
analogue models of gravity. 



3.1 The acoustic geometry in BEC 

Let us start by very briefly reviewing the derivation of the acoustic metric for a BEC 
system, and show that the equations for the phonons of the condensate closely mimic 
the dynamics of a scalar field in a curved spacetime. In the dilute gas approximation, 
one can describe a Bose gas through a quantum field f satisfying 

m^$=(-^-V 2 + V cxt {x) + K(a)&v)$. (1) 

m is the mass of the atoms, a is the scattering length for the atoms and K parametrises 
the strength of the interactions between the different bosons in the gas. It can be re- 
expressed in terms of the scattering length a as 

, . 4nah 2 

K{a) = . (2) 

m 

As usual, the quantum field can be separated into a macroscopic (classical) con- 
densate and a fluctuation: W — y + <p, with (W) = y. Then, by adopting the self- 
consistent mean field approximation 

9 t 99~2(^ t ^) (p + ((p(p) (p\ (3) 

one can arrive at the set of coupled equations: 
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ih—y(t,x) = (-^V 2 + V ext (x) + Kn c ) y(/,3 



-K{2n\if(t,x)+my*(t,x)}; (4) 



ih^- t (p(t,x) = (-^-V 2 + V axt {x) + K2n T )(p(t,x) 



-Km T (p\t,x). (5) 



Here 



n c =|V/(f,x)| 2 ; 


m c = y/ 2 (t,x); 


ft ={$*$); 


m= (<p<p); 


nj = n c + «; 


mr = m c + m. 



(6) 

(7) 
(8) 

In general one will have to solve both equations for y/ and <p simultaneously. The 
equation for the condensate wave function y is closed only when the back-reaction 
effects due to the fluctuations are neglected. (The back-reaction being hidden in the 
quantities m and «.) This approximation leads then to the so-called Gross-Pitaevskii 
equation and can be checked a posteriori to be a good description of dilute Bose- 
Einstein condensates near equilibrium configurations. 

Adopting the Madelung representation for the wave function y of the condensate 



\j/(t,x) = y/n c {t,x) exp[-i6(t,x)/h}, (9) 

and defining an irrotational "velocity field" by v = VO/m, the Gross-Pitaevskii 
equation can be rewritten as a continuity equation plus an Euler equation: 

— n c + V-(n c v)=0, (10) 

4v + vfe + y ext(f ,x) + ,„ c -|l^)Uo. (11) 

at \ 2 2m ^Jn c J 

These equations are completely equivalent to those of an irrotational and inviscid 
fluid apart from the existence of the so-called quantum potential 

Vquantum = -ti 2 V 2 yJTT c / '(2my/n^) , (12) 

which has the dimensions of an energy. 

If we write the mass density of the Madelung fluid as p = m n c , and use the fact 
that the flow is irrotational we can write the Euler equation in the more convenient 
Hamilton-Jacobi form: 

d five} 2 , , Ti 2 V 2 JJ\\ 

at \ 2m 2m ^/n c J 
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When the gradients in the density of the condensate are small one can neglect the 
quantum stress term leading to the standard hydrodynamic approximation. 

Let us now consider the quantum perturbations above the condensate. These can 
be described in several different ways, here we are interested in the "quantum acous- 
tic representation" 



w-^(^*-'f *)■ <i4) 

where n\, 9\ are real quantum fields. By using this representation Equation (0 can 
be rewritten as 

d,ni + -V-(niV0+n c V6 l )=0, (15) 

m \ I 

~ 1 Ti 2 

d,9i + -ve • V0! + «•(«) m - — £> 2 «i = o. (16) 

m 2m 

Here £>2 represents a second-order differential operator obtained from linearizing 
the quantum potential. Explicitly: 

r, - * -3/2 rT7 2/ +l/2y| - . 1 -l/2rf/ -1/2- *, , , _ 

D2ii = --n c [V (n c )J«i + -b c V (n c «!). (17) 

The equations we have just written can be obtained easily by linearizing the Gross- 
Pitaevskii equation around a classical solution: n c — > n c + ni, ^ — > + 0i. It is 
important to realize that in those equations the back-reaction of the quantum fluc- 
tuations on the background solution has been assumed negligible. We also see in 
Equations ([131 [Tol l, that time variations of V ex t and time variations of the scattering 
length a appear to act in very different ways. Whereas the external potential only 
influences the background Equation (TT~3T > (and hence the acoustic metric in the ana- 
logue description), the scattering length directly influences both the perturbation and 
background equations. From the previous equations for the linearised perturbations 
it is possible to derive a wave equation for 9\ (or alternatively, for «i). All we need 
is to substitute in Equation (fT~5T > the n\ obtained from Equation ( fT6] >. This leads to a 
PDE that is second-order in time derivatives but infinite order in space derivatives - 
to simplify things we can construct the symmetric 4x4 matrix 



f» v (t,x) 



f 00 ■■ f j 



r ■ f h 



(18) 



(Greek indices run from 0-3, while Roman indices run from 1-3 .) Then, introducing 
(3+l)-dimensional space-time coordinates 

x* = {t\j) (19) 
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the wave equation for Q\ is easily rewritten as 

W v <3v0i) = O. 



(20) 



Where the f^ v are differential operators acting on space only. Now, if we make a 
spectral decomposition of the field Q\ we can see that for wavelengths larger than 



E, = ft/ mc sound (4 corresponds to the "healing length" and c somi d(a,n c ) 



Jc(rt) n c 



). 



the terms coming from the linearization of the quantum potential (the D2) can be 
neglected in the previous expressions, in which case the f^ v can be approximated 
by scalars, instead of differential operators. Then, by identifying 



UV fjlV 



=F\ 



(21) 



the equation for the field Q\ becomes that of a (massless minimally coupled) quan- 
tum scalar field over a curved background 



AG, 



1 



f v d v )e, = o, 



(22) 



with an effective metric of the form 



^ v (f,x) = 



mc S0UBd (a,n c ) 



-{c S ound(a,«c) 2 -v 2 } 



Si} 



(23) 



Here the magnitude c saan d(n c ,a) represents the speed of the phonons in the medium: 

2 _ K(a) n c 



Csound( fl j w e) 

and v; is the velocity field of the fluid flow, 



v* = -v,e. 

m 



(24) 



(25) 



3.2 Lorentz violation in BEC 



It is interesting to consider the case in which the above "hydrodynamical" approxi- 
mation for BECs does not hold. In order to explore a regime where the contribution 
of the quantum potential cannot be neglected we can use the so called eikonal ap- 
proximation, a high-momentum approximation where the phase fluctuation Q\ is it- 
self treated as a slowly-varying amplitude times a rapidly varying phase. This phase 
will be taken to be the same for both n\ and 8\ fluctuations. In fact, if one discards 
the unphysical possibility that the respective phases differ by a time varying quan- 
tity, any time-independent difference can be safely reabsorbed in the definition of 
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the (complex) amplitudes srfe,srf p . Specifically, we shall write 

0i(f,x)=Re{^ e expH</>)}, 
ni(t,x) = Re{^/ p exp(— /</>)} . 



(26) 
(27) 



As a consequence of our starting assumptions, gradients of the amplitude, and gra- 
dients of the background fields, are systematically ignored relative to gradients of 
0. Note however, that what we are doing here is not quite a "standard" eikonal ap- 
proximation, in the sense that it is not applied directly on the fluctuations of the field 
y/{t,x) but separately on their amplitudes and phases p\ and §\. We can then adopt 
the notation 

(0=-£; ki = Vi<j). (28) 

Then the operator D 2 can be approximated as 



D 2 n\ 



1 



n c k n\. 



(29) 



A similar result holds for D 2 acting on 8\ . That is, under the eikonal approximation 
we effectively replace the operator D 2 by the function 



D 2 ^r — n~ l k 2 . 

2 c 

For the matrix f^ v this effectively results in the replacement 

-l 



f.i _+ 
f iJ -> 



K{a) 

K{a) + 
V'0o 



Ti 2 k 2 
Am n c 

h 2 k 2 

Am n c 

K{a) 



vn 



h 2 k 2 

Am n c 



n c S'J V'flb 



k(cl) 



h 2 k 2 

Am n c 



V-% 



(30) 

(31) 
(32) 
(33) 
(34) 



(As desired, this has the net effect of making f^ v a matrix of numbers, not opera- 
tors.) The physical wave equation (l20t now becomes a nonlinear dispersion relation 



/ uu a 1 + {f +r)co h + f j h kj = o. 



(35) 



After substituting the approximate D 2 into this dispersion relation and rearranging, 
we see (remember: k 2 = ||A:|| 2 = 5' 7 k\ kj) 



co 2 + 2 v' coki - 



n c k 2 



K{a) 



Amn c 



(v i k i ) 2 =0. 



(36) 
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That is (with v^iV^o) 

. 2 n c k 2 



K ( a ) + 1 k 

\mn c 



(CO - Vq kt 
Introducing the speed of sound c SOU nd this takes the form 



(37) 



n 

2m 



<o = v l oki±\c 2 soanii ki+[—ki) . (38) 



2 



We then see that BEC is a paradigmatic framework where a spacetime geometry 
emerges at low energies and Lorentz invariance is as an accidental (never exact) 
symmetry. This symmetry is naturally broken at high energies and appears emi- 
nently in modified dispersion relations for the quasi-particles living above the con- 
densate background. 



4 Modified dispersion relations and their naturalness 

As mentioned before, not only analogue models but also several QG scenarios 
played an important role in motivating searcher for departures from Lorentz invari- 
ance and in most of these models, LV enters through modified dispersion relations 
of the sort (1381 1. These relations can be cast in the general form 

E 2 = p 2 + m 2 +f(E,p;ii;M), (39) 

where the low energy speed of light c = 1 ; E and p are the particle energy and mo- 
mentum, respectively; pi is a particle-physics mass-scale (possibly associated with a 
symmetry breaking/emergence scale) and M denotes the relevant QG scale. Gener- 
ally, it is assumed that M is of order the Planck mass: M ~ M P \ m 1 .22 x 10 19 GeV, 
corresponding to a quantum (or emergent) gravity effect. Note that we assumed a 
preservation of rotational invariance by QG physics and that only boost invariance is 
affected by Planck-scale corrections. This does not need to be the case (see however 
||33ll for a discussion about this assumption) and constraints on possible breakdown 
of rotational invariance have been considered in the literature (especially in the con- 
text of the minimal standard model extension). We assume it here only for simplicity 
and clarity in assessing later the available constraints on the EFT framework. 

Of course, once given (l39l the natural thing to do is to expand the function 
f(E,p;jj,;M) in powers of the momentum (energy and momentum are basically 
indistinguishable at high energies, although they are both taken to be smaller than 
the Planck scale), 

E 2 =/ + m 2 + £%><', (40) 
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where the lowest order LV terms (p, p 2 , /? 3 , p 4 ) have primarily been considered 

an 

About this last point some comments are in order. In fact, from a EFT point of 
view the only relevant operators should be the lowest order ones, i.e. those of mass 
dimension 3,4 corresponding to terms of order p and p 2 in the dispersion relation. 
Situations in which higher order operators "weight" as much as the lowest order 
ones are only possible at the cost of a severe, indeed arbitrary, fine tuning of the 
coefficients fjj. 

However, we do know by now (see further discussion below) that current obser- 
vational constraints are tremendous on dimension 3 operators and very severe on 
dimension 4 ones. This is kind of obvious, given that these operators would end up 
modifying the dispersion relation of elementary particles at low energies. Dimen- 
sion 3 operator would dominate at p — >• while the dimension 4 ones would gener- 
ically induce a, species dependent, constant shift in the limit speed for elementary 
particles. 

Of course one might be content to limit oneself to the study of just these terms 
but we stress that emergent gravity scenarios, e.g. inspired by analogue gravity mod- 
els, or QG gravity models, strongly suggest that if the origin of the breakdown of 
Lorentz invariance is rooted in the UV behaviour of gravitational physics then it 
should be naturally expected to become evident only at high energies. So one would 
then predict a hierarchy of LV coefficients of the sort 

*h = mfp fh = rh^, fc = »fc^, ?)4 = ii4 W _. (4i) 

In characterizing the strength of a constraint one can then refer to the tj„ without 
the tilde, so to compare to what might be expected from Planck-suppressed LV. In 
general one can allow the LV parameters r\i to depend on the particle type, and in- 
deed it turns out that they must sometimes be different but related in certain ways for 
photon polarization states, and for particle and antiparticle states, if the framework 
of effective field theory is adopted. 



4.1 The naturalness problem 

While the above hierarchy (l4TT i might seem now a well motivated framework within 
which asses our investigations, it was soon realized J83l that it is still quite unnatural 
from an EFT point of view. The reason is pretty simple: in EFT radiative corrections 
will generically allow the percolation of higher dimension Lorentz violation to the 
lowest dimension terms due to the coupling and self-couplings of particles ll83l . In 



2 I disregard here the possible appearance of dissipative terms (82) in the dispersion relation, as 
this would correspond to a theoiy with unitarity loss and to a more radical departure from stan- 
dard physics than that envisaged in the framework discussed herein (albeit a priori such dissipative 
scenarios are logically consistent and even plausible within some quantum/emergent gravity frame- 
works). 
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EFT loop integrals will be naturally cut-off at the EFT breaking scale, if such scale 
is as well the Lorentz breaking scale the two will basically cancel leading to un- 
suppressed, couplings dependent, contributions to the propagators. Hence radiative 
corrections will not allow a dispersion relation with only p i or p 4 Lorentz breaking 
terms but will automatically induce extra unsuppressed LV terms in p and p 2 which 
will be naturally dominant. 

Several ideas have been advanced in order to justify such a "naturalness problem" 
(see eg. 11331 ). it would be cumbersome to review here all the proposals, but one can 
clearly see that the most straightforward solution for this problem would consist in 
breaking the degeneracy between the EFT scale and the Lorentz breaking one. This 
can be achieved in two alternative ways. 



4.1.1 A new symmetry 

Most of the aforementioned proposals implicitly assume that the Lorentz breaking 
scale is the Planck scale. One then needs the EFT scale (which can be naively iden- 
tified with what we called previously jj.) to be different from the Planck scale and 
actually sufficiently small so that the lowest order "induced" coefficients can be sup- 
pressed by suitable small rations of the kind ji p /M q where p,q are some positive 
powers. 

A possible solution in this direction can be provided by introducing what is com- 
monly called a "custodial symmetry" something that forbids lower order operators 
and, once broken, suppress them by introducing a new scale. The most plausible 
candidate for this role was soon recognized to be Super Symmetry (SUSY) J84l[85l . 
SUSY is by definition a symmetry relating fermions to bosons i.e. matter with inter- 
action carriers. As a matter of fact, SUSY is intimately related to Lorentz invariance. 
Indeed, it can be shown that the composition of at least two SUSY transformations 
induces space-time translations. However, SUSY can still be an exact symmetry 
even in presence of LV and can actually serve as a custodial symmetry preventing 
certain operators to appear in LV field theories. 

The effect of SUSY on LV is to prevent dimension < 4, renormalizable LV oper- 
ators to be present in the Lagrangian. Moreover, it has been demonstrated 18411851 
that the renormalization group equations for Supersymmetric QED plus the addi- 
tion of dimension 5 LV operators a la Myers & Pospelov 11371 do not generate lower 
dimensional operators, if SUSY is unbroken. However, this is not the case for our 
low energy world, of which SUSY is definitely not a symmetry. 

The effect of soft SUSY breaking was also investigated in [84, 85|. It was found 
there that, as expected, when SUSY is broken the renormalizable operators are gen- 
erated. In particular, dimension k ones arise from the percolation of dimension K + 2 
LV operators^. The effect of SUSY soft-breaking is, however, to introduce a sup- 
pression of order m 2 /Mp\ (k = 3) or (m s /Mpi) 2 (k = 4), where m s ~ 1 TeV is the 
scale of SUSY soft breaking. Although, given present constraints, the theory with 



We consider here only K" = 3 , 4, for which these relationships have been demonstrated. 
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K = 3 needs a lot of fine tuning to be viable, since the SUSY-breaking-induced sup- 
pression is not enough powerful to kill linear modifications in the dispersion relation 
of electrons, if K = 4 then the induced dimension 4 terms are suppressed enough, 
provided m s < 100 TeV. Current lower bounds from the Large Hadron Collider are 
at most around 950 GeV for the most simple models of SUSY ll86l (the so called 
"constrained minimal supersymmetric standard model", CMSSM). 

Finally, it is also interesting to note that the analogue model of gravity can be used 
as a particular implementation of the above mentioned mechanism for avoiding the 
so called naturalness problem via a custodial symmetry. This was indeed the case of 
multi-BEC (831881. 



4.1.2 Gravitational confinement of Lorentz violation 

The alternative to the aforementioned scenario is to turn the problem upside down. 
One can in fact assume that the Lorentz breaking scale (the M appearing in the 
above dispersion relations) is not set by the Planck scale while the latter is the EFT 
breaking scale. If in addition one starts with a theory which has higher order Lorentz 
violating operators only in the gravitational sector, then one can hope that the grav- 
itational coupling Gn ~ Mpj 2 will let them "percolate" to the matter sector however 
it will do so introducing factors of the order (M/Mp\) 2 which can become strong 
suppression factors if M -C Mp\. This is basically the idea at the base of the work 
presented in J89l which applies it to the special case of Hofava-Lifshitz gravity. 
There it was shown that indeed a workable low energy limit of the theory can be 
derived through this mechanism which apparently is fully compatible with extant 
constraints on Lorentz breaking operators in the matter sector. We think that this 
new route deserves further attention and should be more deeply explored in the fu- 
ture. 



5 Dynamical frameworks 

Missing a definitive conclusion about the naturalness problem, the study of LV 
theories has basically proceeded by considering separately extensions of the Stan- 
dard Model based on naively power counting renormalizable operators or non- 
renormalizable operators (at some given mass dimension). In what follows we shall 
succinctly describe these frameworks before to discuss theoretical alternatives. 



5.1 SME with renormalizable operators 

Most of the research in EFT with only renormalizable (i.e. mass dimension 3 and 4) 
LV operators has been carried out within the so called (minimal) SME l22l . It con- 
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sists of the standard model of particle physics plus all Lorentz violating renormal- 
izable operators (i.e. of mass dimension < 4) that can be written without changing 
the field content or violating gauge symmetry. The operators appearing in the SME 
can be conveniently classified according to their behaviour under CPT. Since the 
most common particles used to cast constraints on LV are photons and electrons, a 
prominent role is played by LV QED. 

If we label by ± the two photon helicities, we can write the photon dispersion 
relation as Ir90l 

E = (l+p±o)\p\ (42) 

where p and a depend on LV parameters appearing in the LV QED Lagrangian, 
as defined in J7). Note that the dependence of the dispersion relation on the photon 
helicity is due to the fact that the SME generically also contemplates the possibility 
of a breakdown of rotational invariance. 

We already gave (see section |4]i motivations for assuming rotation invariance 
to be preserved, at least in first approximation, in LV contexts. If we make this 
assumption, we obtain a major simplification of our framework, because in this case 
all LV tensors must reduce to suitable products of a time-like vector field, which 
is usually called u a and, in the preferred frame, is assumed to have components 
(1,0,0,0). Then, the rotational invariant LV operators are 

1 _ ° 1 ° 

-^« M V / 757 M V / + TjicitfiUvWY^ D v \j/+ -idunU v \j/Y5r ® V V (43) 



for electrons and 



l -{k F )u K r] Xil u v F KX F^ (44) 



for photons. 

The high energy {Mp\ ^> E ^> m) dispersion relations for QED can be expressed 
as (see Q and references therein for more details) 

E\ = m] + p 2 +fPp +f^ 2) p 2 electrons (45) 

E y = ( 1 + /y 2) )P 2 Photons (46) 

where f e = —2bs, f e ' = — (c — ds), and fy ' = Icf /2 with s = ±1 the helicity state 
of the electron |7|. The positron dispersion relation is the same as d45l l replacing 
p —¥ —p, this will change only the / e term. 

We notice here that the typical energy at which a new phenomenology should 

(2) 

start to appear is quite low. In fact, taking for example fe ~ 0(1), one finds 
that the corresponding extra-term is comparable to the electron mass m precisely at 
p ~ m ~ 5 1 1 keV. Even worse, for the linear modification to the dispersion relation, 
we would have, in the case in which /j ~ 0(1), that p t h ~ m 2 /Mp\ ~ 10~ 17 eV. 
(Notice that this energy corresponds by chance to the present upper limit on the 
photon mass, niy < 10~ 18 eV ||9T1 .) As said, this implied strong constraints on the 
parameters and was a further mortivation for exploring the QG preferred possibil- 
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ity of higher order Lorentz violating operators and consequently try to address the 
naturalness problem. 



5.2 Dimension five operators SME 

An alternative approach within EFT is to study non-renormalizable operators. 
Nowadays it is widely accepted that the SM could just be an effective field theory 
and in this sense its renormalizability is seen as a consequence of neglecting some 
higher order operators which are suppressed by some appropriate mass scale. It is a 
short deviation from orthodoxy to imagine that such non-renormalizable operators 
can be generated by quantum gravity effects (and hence be naturally suppressed 
by the Planck mass) and possibly associated to the violation of some fundamental 
space-time symmetry like local Lorentz invariance. 

Myers & Pospelov J37l found that there are essentially only three operators of 
dimension five, quadratic in the fields, that can be added to the QED Lagrangian 
preserving rotation and gauge invariance, but breaking local Lo 

These extra-terms, which result in a contribution of 0(E/Mp\) to the dispersion 
relation of the particles, are the following: 

P -u m F ma (u-d){u n F n ") + -^-u m W7 m {^ + ^75){u-d) 2 W, (47) 



2M P1 '"" v " " ' 2M P . 

where F is the dual of F and |, £12 are dimensionless parameters. All these terms 
also violate the CPT symmetry. More recently, this construction has been extended 
to the whole SM <92l. 

From d47T > the dispersion relations of the fields are modified as follows. For the 
photon one has 

(ol=k 2 ±^-k\ (48) 

(the + and — signs denote right and left circular polarisation), while for the fermion 
(with the + and — signs now denoting positive and negative helicity states) 

Ei=p 2 + m 2 + r, ± {-, (49) 

with rj-t = 2(£i ± £2). For the antifermion, it can be shown by simple "hole in- 
terpretation" arguments that the same dispersion relation holds, with 77^ = — 77^ 
where af and / superscripts denote respectively anti-fermion and fermion coeffi- 
cients 



4 Actually these criteria allow the addition of other (CPT even) terms, but these would not lead 
to modified dispersion relations (they can be thought of as extra, Planck suppressed, interaction 
terms) (92). 
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As we shall see, observations involving very high energies can thus potentially 
cast 0(1) and stronger constraint on the coefficients defined above. A natural ques- 
tion arises then: what is the theoretically expected value of the LV coefficients in the 
modified dispersion relations shown above? 

This question is clearly intimately related to the meaning of any constraint pro- 
cedure. Indeed, let us suppose that, for some reason we do not know, because we do 
not know the ultimate high energy theory, the dimensionless coefficients 77'"', that 
in principle, according to the Dirac criterion, should be of order (9(1), are defined 
up to a dimensionless factor of m e /Mp\ ~ 10~ 22 . (This could well be as a result of 
the integration of high energy degrees of freedom.) Then, any constraint of order 
larger than 10~ 22 would be ineffective, if our aim is learning something about the 
underlying QG theory. 

This problem could be further exacerbated by renormalization group effects, 
which could, in principle, strongly suppress the low-energy values of the LV co- 
efficients even if they are 0( 1 ) at high energies. Let us, therefore, consider the evo- 
lution of the LV parameters first. Bolokhov & Pospelov J92l addressed the problem 
of calculating the renormalization group equations for QED and the Standard Model 
extended with dimension-five operators that violate Lorentz Symmetry. 

In the framework defined above, assuming that no extra physics enters between 
the low energies at which we have modified dispersion relations and the Planck 
scale at which the full theory is defined, the evolution equations for the LV terms in 
Eq. d4Tb that produce modifications in the dispersion relations, can be inferred as 

dC,\ 25 a dC, 2 25 a 5 a g dt, \ a 2a 

dt 12 % dt 12 % 12 % dt 12 % 3 % 

where a = e 2 j\n ~ 1/137 (h = 1) is the fine structure constant and t = ln(^ 2 //Xg) 
with /i and /Xrj two given energy scales. (Note that the above formulae are given 
to lowest order in powers of the electric charge, which allows one to neglect the 
running of the fine structure constant.) 

These equations show that the running is only logarithmic and therefore low 
energy constraints are robust: 0(1) parameters at the Planck scale are still 0(1) at 
lower energy. Moreover, they also show that 77+ and T7_ cannot, in general, be equal 
at all scales. 



5.3 Dimension six operators SME 

If CPT ends up being a fundamental symmetry of nature it would forbid all of the 
above mentioned operators (hence pushing at further high energies the emergence 
of Lorentz breaking physics). It makes then sense to consider dimension six, CPT 
even, operators which furthermore do give rise to dispersion relations of the kind 
appearing in the above mentioned BEC analogue gravity example. 
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The CPT even dimension 6 LV terms have only recently been computed 
through the same procedure used by Myers & Pospelov for dimension 5 LV. The 
known fermion operators are 

-^W( u -D) 2 (a { L 5) P L + a { R 5) P R ) W 
-y?(«-0) 3 (" • Y){af ] P L + a R 6) P R )y (51) 



A/p, 



T xj7{u-D)D(u ■ y)(a { L 6) P L + a R 6) P R )w , 



where Pr± are the usual left and right spin projectors Pr.l = (1 ± y 5 )/2 and where D 
is the usual QED covariant derivative. All coefficients a are dimensionless because 
we factorize suitable powers of the Planck mass. 
The known photon operator is 

- ^P^F^u^iu ■ d)F GV . (52) 

From these operators, the dispersion relations of electrons and photons can be 
computed, yielding 



*-/^ = ^rf + 4V + 4V|r 



a (6) £ 3 a (6) £ 3 

+^^( £ + ^) + -^( £ -^) (53) 

(0 2_ k 2 =j3 (6)]S 

where m is the electron mass and where s = a • p/|p| is the helicity of the electrons. 
Also, notice that a term proportional to E is generated. 

Because the high-energy fermion states are almost exactly chiral, we can further 
simplify the fermion dispersion relation eq. (l54t (we pose R = +, L = — ) 

£ 2 =„ 2 + , w 2 +/j 4 y+/l 6) ^. (55) 

Being suppressed by a factor of order m/Mp\, we will drop in the following the 
quadratic contribution f± p 2 , indeed this can be safely neglected, provided that 
E > ^/mMpi . Let me stress however, that this is exactly an example of a dimension 4 
LV term with a natural suppression, which for electron is of order m e /Mp\ ~ 10 -22 . 
Therefore, any limit larger than 10~ 22 placed on this term would not have to be 
considered as an effective constraint. To date, the best constraint for a rotational 
invariant electron LV term of dimension 4 is <9(10~ 16 ) [94|. 

Coming back to equation (I5"5l l. it may seem puzzling that in a CPT invariant 
theory we distinguish between different fermion helicities. However, although they 
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are CPT invariant, some of the LV terms displayed in Eq. ((52) are odd under P and 
T. 

CPT invariance allows us to determine a relationship between the LV coefficients 
of the electrons and those of the positrons. Indeed, to obtain these we must consider 
that, by CPT, the dispersion relation of the positron is given by d54l i. with the re- 
placements s — >• — s and p —> —p. This implies that the relevant positron coefficients 
/oositron are sucn ^ a * / + = / - > where e± indicates a positron of positive/negative 
helicity (and similarly for the e±). 



5.4 Other frameworks 

Picking up a well defined dynamical framework is sometimes crucial in discussing 
the phenomenology of Lorentz violations. In fact, not all the above mentioned "win- 
dows on quantum gravity" can be exploited without adding additional information 
about the dynamical framework one works with. Although cumulative effects ex- 
clusively use the form of the modified dispersion relations, all the other "windows" 
depend on the underlying dynamics of interacting particles and on whether or not 
the standard energy-momentum conservation holds. Thus, to cast most of the con- 
straints on dispersion relations of the form (140) , one needs to adopt a specific theo- 
retical framework justifying the use of such deformed dispersion relations. 

The previous discussion mainly focuses on considerations based on Lorentz 
breaking EFTs. This is indeed a conservative framework within which much can be 
said (e.g. reaction rates can still be calculated) and from an analogue gravity point of 
view it is just the natural frame to work within. Nonetheless, this is of course not the 
only dynamical framework within which a Lorentz breaking kinematics can be cast. 
Because the EFT approach is nothing more than a highly reasonable, but rather ar- 
bitrary "assumption", it is worth studying and constraining additional models, given 
that they may evade the majority of the constraints discussed in this review. 



5.4.1 D-brane models 

We consider here the model presented in 115211541 . in which modified dispersion rela- 
tions are found based on the Liouville string approach to quantum space-time J95l . 
Liouville-string models of space-time foam |{95l motivate corrections to the usual 
relativistic dispersion relations that are first order in the particle energies and that 
correspond to a vacuum refractive index 77 = 1 — (E/Mp\) a , where a = 1. Models 
with quadratic dependences of the vacuum refractive index on energy: a = 2 have 
also been considered {59). 

In particular, the D-particle realization of the Liouville string approach predicts 
that only gauge bosons such as photons, not charged matter particles such as elec- 
trons, might have QG-modified dispersion relations. This difference may be traced 
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to the facts that 11961 excitations which are charged under the gauge group are rep- 
resented by open strings with their ends attached to the D-brane ll97l , and that only 
neutral excitations are allowed to propagate in the bulk space transverse to the brane. 
Thus, if we consider photons and electrons, in this model the parameter 77 is forced 
to be null, whereas £, is free to vary. Even more importantly, the theory is CPT even, 
implying that vacuum is not birefringent for photons (§+ = <*_). 



5.4.2 Doubly Special Relativity 

Lorentz invariance of physical laws relies on only few assumptions: the principle 
of relativity, stating the equivalence of physical laws for non-accelerated observers, 
isotropy (no preferred direction) and homogeneity (no preferred location) of space- 
time, and a notion of precausality, requiring that the time ordering of co-local events 
in one reference frame be preserved M\ W\ USD EECI HQ2] [[03] QUI HSU ■ 

All the realizations of LV we have discussed so far explicitly violate the principle 
of relativity by introducing a preferred reference frame. This may seem a high price 
to pay to include QG effects in low energy physics. For this reason, it is worth 
exploring an alternative possibility that keeps the relativity principle but that relaxes 
one or more of the above postulates. 

For example, relaxing the space isotropy postulate leads to the so-called Very 
Special Relativity framework 111061 . which was later on understood to be described 
by a Finslerian-type geometry H1071 11081 11091 . In this example, however, the gen- 
erators of the new relativity group number fewer than the usual ten associated with 
Poincare invariance. Specifically, there is an explicit breaking of the 0(3) group 
associated with rotational invariance. 

One may wonder whether there exist alternative relativity groups with the same 
number of generators as special relativity. Currently, we know of no such generaliza- 
tion in coordinate space. However, it has been suggested that, at least in momentum 
space, such a generalization is possible, and it was termed "doubly" or "deformed" 
(to stress the fact that it still has 10 generators) special relativity, DSR. Even though 
DSR aims at consistently including dynamics, a complete formulation capable of 
doing so is still missing, and present attempts face major problems. Thus, at present 
DSR is only a kinematic theory. Nevertheless, it is attractive because it does not 
postulate the existence of a preferred frame, but rather deforms the usual concept of 
Lorentz invariance in the following sense. 

Consider the Lorentz algebra of the generators of rotations , L,, and boosts, Bf. 

[Li,Lj] = 1 e t jk L k ; [Li,Bj] = 1 £, ;/t B k ; [Bi,Bj] = — ; e t j k L k (56) 

(Latin indices i,j, . . . run from 1 to 3) and supplement it with the following com- 
mutators between the Lorentz generators and those of translations in spacetime (the 
momentum operators Pq and Pf): 

[L„Po]=0; [Li,Pj] = ie ijk P k ; (57) 
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[B i ,P ]=ifi(^JP i ; (58) 



[Bi,Pj] 



P\„ , r fP\ PjPj 



8ijfa[-)P + fa 



K \K K 



(59) 



where K is some unknown energy scale. Finally, assume [Pi,Pj] = 0. The commuta- 
tion relations J58l-(l59l are given in terms of three unspecified, dimensionless struc- 
ture functions fa, fa, and fa, and they are sufficiently general to include all known 
DSR proposals — the DSR1 OH, DSR2 0TUQ31, and DSR3 EB). Further- 
more, in all the DSRs considered to date, the dimensionless arguments of these 
functions are specialized to 



so rotational symmetry is completely unaffected. For the k —$■ +°° limit to reduce to 
ordinary special relativity, fa and fa must tend to 1, and fa must tend to some finite 
value. 

DSR theory postulates that the Lorentz group still generates space-time symme- 
tries but that it acts in a non-linear way on the fields, such that not only is the speed 
of light c an invariant quantity, but also that there is a new invariant momentum scale 
K which is usually taken to be of the order of Mpi. Note that DSR-like features are 
found in models of non-commutative geometry, in particular in the iC-Minkowski 
framework 0114111151 . as well as in non-canonically non commutative field theories 

ma. 

Concerning phenomenology, an important point about DSR in momentum space 
is that in all three of its formulations (DSR1 fljij], DSR2 EUEE), and DSR3 EH) 
the component of the four momentum having deformed commutation with the boost 
generator can always be rewritten as a non-linear combination of some energy- 
momentum vector that transforms linearly under the Lorentz group 111 171 . For ex- 
ample in the case of DSR2 11111111211 one can write s 

(61) 



I-TTo/K ' 

Pi = : - , ■ (62) 

1-flb/K 

It is easy to ensure that while n satisfies the usual dispersion relation 7Iq — tz 2 = m 2 
(for a particle with mass m), E and pi satisfy the modified relation 

( 1 - to 2 / K 2 ) E 2 + 2 K~ ' m 2 E-p 2 = m 2 . (63) 

Furthermore, a different composition for energy-momentum now holds, given that 
the composition for the physical DSR momentum p must be derived from the stan- 
dard energy-momentum conservation of the pseudo-variable n and in general im- 
plies non-linear terms. A crucial point is that due to the above structure if a threshold 
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reaction is forbidden in relativistic physics then it is going to be still forbidden by 
DSR. Hence many constraints that apply to EFT do not apply to DSR. 

Despite its conceptual appeal, DSR is riddled with many open problems. First, if 
DSR is formulated as described above — that is, only in momentum space — then 
it is an incomplete theory. Moreover, because it is always possible to introduce the 
new variables 7tyi, on which the Lorentz group acts in a linear manner, the only way 
that DSR can avoid triviality is if there is some physical way of distinguishing the 
pseudo-energy e = — 7Cq from the true-energy E, and the pseudo-momentum n from 
the true-momentum p. If not, DSR is no more than a nonlinear choice of coordinates 
in momentum space. 

In view of the standard relations E o ihd t and p -H> —ihV (which are presumably 
modified in DSR), it is clear that to physically distinguish the pseudo-energy e from 
the true-energy E, and the pseudo-momentum n from the true-momentum p, one 
must know how to relate momenta to position. At a minimum, one needs to develop 
a notion of DSR spacetime. 

In this endeavor, there have been two distinct lines of approach, one presuming 
commutative spacetime coordinates, the other attempting to relate the DSR feature 
in momentum space to a non commutative position space. In both cases, several 
authors have pointed out major problems. In the case of commutative spacetime co- 
ordinates, some analyses have led authors to question the triviality [118] or internal 
consistency 111191 11201 11211 of DSR. On the other hand, non-commutative propos- 
als l57l are not yet well understood, although intense research in this direction is 



under way 11221 . Finally we cannot omit the recent development of what one could 
perhaps consider a spin-off of DSR that is Relative Locality, which is based on the 
idea that the invariant arena for classical physics is a curved momentum space rather 
than spacetime (the latter being a derived concept) 01231 . 

DSR and Relative Locality are still a subject of active research and debate (see 
e.g. fl 1241 1 1251 1 1261 [1271 ): nonetheless, they have not yet attained the level of ma- 
turity needed to cast robust constraints^. For these reasons, in the next sections we 
focus upon LV EFT and discuss the constraints within this framework. 



6 Experimental probes of low energy LV: Earth based 
experiments 

The world as we see it seems ruled by Lorentz invariance to a very high degree. 
Hence, when seeking tests of Lorentz violations one is confronted with the chal- 
lenge to find or very high precision experiments able to test Special Relativity or 
observe effects which might be sensitive to tiny deviations from standard LI. Within 



5 Note however, that some knowledge of DSR phenomenology can be obtained by considering 
that, as in Special Relativity, any phenomenon that implies the existence of a preferred reference 
frame is forbidden. Thus, the detection of such a phenomenon would imply the falsification of 
both special and doubly-special relativity. An example of such a process is the decay of a massless 
particle. 
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the ansatz we lied down in the previous sections it is clear that the first route is 
practical only when dealing with low energy violations of Lorentz invariance as 
systematically described by the minimal Standard Model extension (mSME) while 
astrophysical tests, albeit much less precise, are the choice too for testing LV in- 
duced by higher order operators. Let us then briefly review the main experimental 
tools used so far in order to perform precision tests of Lorentz invariance in labora- 
tory (for more details see e.g. IT7] l24][T28l ). 



6.1 Penning traps 

In a Penning trap a charged particle can be localized for long times using a com- 
bination of static magnetic and electric fields. Lorentz violating tests are based on 
monitoring the particle cyclotron motion in the magnetic field and Larmor preces- 
sion due to the spin. In fact the relevant frequencies for both these motions are 
modified in the mSME and Penning traps can be made very sensitive to differences 
in these frequencies. 



6.2 Clock comparison Experiments 

Clock comparison experiments are generally performed by considering two atomic 
transition frequencies (which can be considered as two clocks) in the same point 
in space. The basic idea is that as the clocks move in space, they pick out different 
components of the Lorentz violating tensors in the mSME. This would yield a side- 
real drift between the two clocks. Measuring the difference between the frequencies 
over long periods, allows to cast very high precision limits on the parameters in the 
mSME (generally for for protons and neutrons.) 



6.3 Cavity Experiments 

In cavity experiments one casts constraints on the variation of the cavity resonance 
frequency as its orientation changes in space. While this is intrinsically similar to 
clock comparison experiments, these kind of experiments allows to cast constraints 
also on the electromagnetic sector of the mSME. 
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6.4 Spin polarized torsion balance 

The electron sector of the mSME can be effective constrained via spin-torsion bal- 
ances. An example is an octagonal pattern of magnets which is constructed so to 
have an overall spin polarization in the octagons plane. Four of these octagons 
are suspended from a torsion fiber in a vacuum chamber. This arrangement of the 
magnets give an estimated net spin polarization equivalent to « 10 23 aligned elec- 
tron spins. The whole apparatus is mounted on a turntable. As the turntable moves 
Lorentz violation in the mSME produces an interaction potential for non-relativistic 
electrons which induces a torque on the torsion balance. The torsion fiber is then 
twisted by an amount related to the relevant LV coefficients. 



6.5 Neutral mesons 

In the mSME one expects an orientation dependent change in the mass difference 
e.g. of neutral Kaons. By looking for sidereal variations or other orientation effects 
one can derive bounds on each component of the relevant LV coefficients. 



7 Observational probes of high energy LV: astrophysical QED 
reactions 

Let us begin with a brief review of the most common types of reaction exploited in 
order to give constraints on the QED sector. 

For definiteness, we refer to the following modified dispersion relations: 

E; = k 2 + 1, { " ] —^ Photon (64) 

E 2 , = m 2 e + p 2 + 77 i n) -^-y Electron-Positron , (65) 

where, in the EFT case, we have §M = £J n) = (-)"|i n) and 77 W = Tji B) = (-)"T7 W . 



7.1 Photon time of flight 

Although photon time-of-flight constraints currently provide limits several orders of 
magnitude weaker than the best ones, they have been widely adopted in the astro- 
physical community. Furthermore they were the first to be proposed in the seminal 
paper 11251 . More importantly, given their purely kinematic al nature, they may be 
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applied to a broad class of frameworks beyond EFT with LV. For this reason, we 
provide a general description of time-of-flight effects, elaborating on their applica- 
tion to the case of EFT below. 

In general, a photon dispersion relation in the form of (l64t implies that photons 
of different colors (wave vectors k\ and £2) travel at slightly different speeds. Let 
us first assume that there are no birefringent effects, so that <i;| = t, . Then, upon 
propagation on a cosmological distance d, the effect of energy dependence of the 
photon group velocity produces a time delay 

1 l*n— 2 1JI— 2 

*M = !!zi*L^_ § W rf| (66 ) 

which clearly increases with d and with the energy difference as long as n > 2. 
The largest systematic error affecting this method is the uncertainty about whether 
photons of different energy are produced simultaneously in the source. 

So far, the most robust constraints on B, ( 3 ', derived from time of flight differences, 
have been obtained within the D— brane model (discussed in section IB". 4. jj from a 
statistical analysis applied to the arrival times of sharp features in the intensity at 
different energies from a large sample of GRBs with known redshifts 01291 , leading 
to limits tip 3 ' < <9(10 3 ). A recent example illustrating the importance of systematic 
uncertainties can be found in [ 130|, where the strongest limit £, ( 3 ' < 47 is found by 
looking at a very strong flare in the TeV band of the AGN Markarian 501. 

One way to alleviate systematic uncertainties — available only in the context 
of birefringent theories, such as the one with n = 3 in EFT — would be to mea- 
sure the velocity difference between the two polarization states at a single energy, 
corresponding to 

At = 2\^ ) \kd/M n . (67) 

This bound would require that both polarizations be observed and that no spurious 
helicity-dependent mechanism (such as, for example, propagation through a bire- 
fringent medium) affects the relative propagation of the two polarization states. 

Let us stress that Eq. d66l l is no longer valid in birefringent theories. In fact, pho- 
ton beams generally are not circularly polarized; thus, they are a superposition of 
fast and slow modes. Therefore, the net effect of this superposition may partially or 
completely erase the time-delay effect. To compute this effect on a generic photon 
beam in a birefringent theory, let us describe a beam of light by means of the as- 
sociated electric field, and let us assume that this beam has been generated with a 
Gaussian width 

E = A ( e '(^r-k + (n Q )z) e -(z-vp) 2 sni i+ + e i(n 0t -k-(n )z) e -(z-vp) 2 sn2 i _\ ^ (6g) 



where Qo is the wave frequency, 8Q0 is the gaussian width of the wave, k (A)) 
is the "momentum" corresponding to the given frequency according to d64l i and 
e± = {e\ ± ie2)/V2 are the helicity eigenstates. Note that by complex conjugation 
e* + =e_. Also, note that ^(a)) = CO^^CO 2 /Mpi. Thus, 
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E=Ae in (t-z) f e ^ni/M Pl z e -(z-vp) 2 snl )i+ + e -^ni/M m z e -(z-vp) 2 sni i \ (69) 
The intensity of the wave beam can be computed as 

E-E* = \A\ 2 f e 2^/M Pl z + e -2i^ /M Pl z\ e -5i2 2 ((z-v+0 2 +(z-v-r) 2 ) 



= 2|A|V 2 ^o(-0 2 cos 2^Qoz) e- 2 ^ {Sn °' )2 . (70) 

V m p\ J 

This shows that there is an effect even on a linearly-polarised beam. The effect is 
a modulation of the wave intensity that depends quadratically on the energy and 
linearly on the distance of propagation. In addition, for a gaussian wave packet, 
there is a shift of the packet centre, that is controlled by the square of B, "'/Mpi and 
hence is strongly suppressed with respect to the cosinusoidal modulation. 



7.2 Vacuum Birefringence 

The fact that electromagnetic waves with opposite "helicities" have slightly different 
group velocities, in EFT LV with n = 3, implies that the polarisation vector of a 
linearly polarised plane wave with energy k rotates, during the wave propagation 
over a distance d, through the angle ll33l Lj 

«M*)-M*) , gg( 3)*rf (72) 

L ZMp] 

Observations of polarized light from a distant source can then lead to a constraint 
on |<i;( 3 '| that, depending on the amount of available information — both on the 
observational and on the theoretical (i.e. astrophysical source modeling) side — can 
be cast in two different ways 11311 : 

1 . Because detectors have a finite energy bandwidth, Eq. (l72t is never probed in real 
situations. Rather, if some net amount of polarization is measured in the band 
k\ < E < k%, an order-of-magnitude constraint arises from the fact that if the 
angle of polarization rotation (T72T > differed by more than n/2 over this band, the 
detected polarization would fluctuate sufficiently for the net signal polarization 
to be suppressed 013211931 , From (l72l , this constraint is 



6 Note that for an object located at cosmological distance (let z be its redshift), the distance d 
becomes 

d(z) = -z\- [' , - dz' , (71) 

#o Jo ^/n A +Q,„(l+zy 

where d(z) is not exactly the distance of the object as it includes a (1 +z) 2 factor in the integrand 
to take into account the redshift acting on the photon energies. 
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£(3)< K -MH (73) 

This constraint requires that any intrinsic polarization (at source) not be com- 
pletely washed out during signal propagation. It thus relies on the mere detection 
of a polarized signal; there is no need to consider the observed polarization de- 
gree. A more refined limit can be obtained by calculating the maximum observ- 
able polarization degree, given the maximum intrinsic value fl 1 331 : 



7T(^)= H(0)^(cos(2e))^ + (sin(20)) 2 ^, (74) 

where 17(0) is the maximum intrinsic degree of polarization, is defined in 
Eq. (l72l and the average is weighted over the source spectrum and instrumental 
efficiency, represented by the normalized weight function &(k) 11321 . Conser- 
vatively, one can set n(0) = 100%, but a lower value may be justified on the 
basis of source modeling. Using (|74l , one can then cast a constraint by requiring 
n(£, ) to exceed the observed value. 
2. Suppose that polarized light measured in a certain energy band has a position an- 
gle o bs with respect to a fixed direction. At fixed energy, the polarization vector 
rotates by the angle (l72ll 7 l: if the position angle is measured by averaging over a 
certain energy range, the final net rotation (A 9) is given by the superposition of 
the polarization vectors of all the photons in that range: 

M2^ = ^, (75) 

(cos (20))^, 

where is given by (l72l . If the position angle at emission 0; in the same energy 
band is known from a model of the emitting source, a constraint can be set by 
imposing 

tan(2(40))<tan(20 obs -20i). (76) 

Although this limit is tighter than those based on eqs. (1731 and (l74l , it clearly 
hinges on assumptions about the nature of the source, which may introduce sig- 
nificant uncertainties. 

In conclusion the fact that polarised photon beams are indeed observed from 
distant objects imposes strong constraints on LV in the photon sector (i.e. on <^' 3 '), 
as we shall see later on. 



7.3 Threshold reactions 

An interesting phenomenology of threshold reactions is introduced by LV in EFT; 
also, threshold theorems can be generalized 11321 . Sticking to the present case of 



7 Faraday rotation is negligible at high energies. 
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rotational invariance and monotonic dispersion relations (see II 1341 for a general- 
ization to more complex situations), the main conclusions of the investigation into 
threshold reactions are that ||3T| 

• Threshold configurations still corresponds to head-on incoming particles and par- 
allel outgoing ones 

• The threshold energy of existing threshold reactions can shift, and upper thresh- 
olds (i.e. maximal incoming momenta at which the reaction can happen in any 
configuration) can appear 

• Pair production can occur with unequal outgoing momenta 

• New, normally forbidden reactions can be viable 

LV corrections are surprisingly important in threshold reactions because the LV 
term (which as a first approximation can be considered as an additional mass term) 
should be compared not to the momentum of the involved particles, but rather to the 
(invariant) mass of the particles produced in the final state. Thus, an estimate for the 
threshold energy is 

fm 2 M'- 2 V /n 
P ^{-^) ' (7?) 

where m is the typical mass of particles involved in the reaction. Interesting values 
for /? t h are discussed, e.g., in |[3T1 and given in Tab.Q] Reactions involving neutrinos 



Table 1 Values of p t h> according to eq. 1771 . for different particles involved in the reaction: neu- 
trinos, electrons and proton. Here we assume T] '"' ~ 1 . 





m v ~0.1 eV 


m e ~ 0.5 MeV 


m p ~ 1 GeV 


n = 2 


0.1 eV 


0.5 MeV 


1 GeV 


n = 3 


500 MeV 


14 TeV 


2PeV 


n = 4 


33 TeV 


74PeV 


3EeV 



are the best candidate for observation of LV effects, whereas electrons and positrons 
can provide results for n = 3 theories but can hardly be accelerated by astrophysical 
objects up to the required energy for n = 4. In this case reactions of protons can be 
very effective, because cosmic -rays can have energies well above 3 EeV. Let us now 
briefly review the main reaction used so far in order to casts constraints. 



7.3.1 LV-allowed threshold reactions: y-decay 



The decay of a photon into an electron/positron pair is made possible by LV because 
energy-momentum conservation may now allow reactions described by the basic 
QED vertex. This process has a threshold that, if £, ~ and n = 3, is set by the 
condition l33l 

(78) 



*tt = (6^ll^lf/|7l£ 3) |) 1 /3. 
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Noticeably, as already mentioned above, the electron-positron pair can now be cre- 
ated with slightly different outgoing momenta (asymmetric pair production). Fur- 
thermore, the decay rate is extremely fast above threshold [1331 and is of the order of 
(lOns)" 1 (n = 3) or (10~ 6 ns)" 1 (n = 4). 



7.3.2 LV-allowed threshold reactions: Vacuum Cerenkov and Helicity Decay 

In the presence of LV, the process of Vacuum Cerenkov (VC) radiation e — >• e ± y 
can occur. If we set B, ~ and n = 3, the threshold energy is given by 

PVC = (m 2 e M/2ri {3) ) 1/3 ~ 11 TeV TJ -1 / 3 . (79) 

Just above threshold this process is extremely efficient, with a time scale of order 
Tvc ~ 10- 9 s J33]. 

A slightly different version of this process is the Helicity Decay (HD, e T — > e±y). 
If TJ+ ^ r\ , an electron/positron can flip its helicity by emitting a suitably polarized 
photon. This reaction does not have a real threshold, but rather an effective one 1331 
— /?hd = (m%M/A Tj) 1 / 3 , where Ai] = \r]\. — r]_ '\ — at which the decay lifetime 
Xhd is minimized. For Ar\ rs 0(1) this effective threshold is around 10 TeV. Note 
that below threshold t H d > Arj~ 3 (p/10 TeV) -8 10~ 9 s, while above threshold t H d 
becomes independent of A r\ |33l . 

Apart from the above mentioned examples of reactions normally forbidden and 
now allowed by LV dispersion relations, one can also look for modifications of nor- 
mally allowed threshold reactions especially relevant in high energy astrophysics. 



7.3.3 LV-allowed threshold reactions: photon splitting and lepton pair 
production 

It is rather obvious that once photon decay and vacuum Cerenkov are allowed also 
the related relations in which respectively the our going lepton pair is replaced by 
two or more photons, y — > 2 7 and 7 — > 37, etc. , or the outgoing photons is replaced 
by an electron-positron pair, e~ — > e~e~e + , are also allowed. 



Photon splitting 

This is forbidden for |W < while it is always allowed if |W > El- When 
allowed, the relevance of this process is simply related to its rate. The most relevant 
cases are 7—)- 77 and 7— >• 37, because processes with more photons in the final state 
are suppressed by more powers of the fine structure constant. 

The 7 — >• 77 process is forbidden in QED because of kinematics and C-parity 
conservation. In LV EFT neither condition holds. However, we can argue that this 
process is suppressed by an additional power of the Planck mass, with respect to 
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y — > 3 7. In fact, in LI QED the matrix element is zero due to the exact cancellation 
of fermionic and anti-fermionic loops. In LV EFT this cancellation is not exact and 
the matrix element is expected to be proportional to at least (%E/Mp{) p , p > 0, as it 
is induced by LV and must vanish in the limit M?\ —> °°. 

Therefore we have to deal only with y— s- 3y. This process has been studied in 
||3T1|1351 . In particular, in 111351 it was found that, if the "effective photon mass" 
m 2 , = t^Ey/Mp^ <C m 2 ., then the splitting lifetime of a photon is approximately 

T n=3 ~ 0.025 4 ~ 5 /~ ' (50 TeV/Ey) s, where / is a phase space factor of order 1 . 
This rate was rather higher than the one obtained via dimensional analysis in IBTI 
because, due to integration of loop factors, additional dimensionless contributions 
proportional to m\ enhance the splitting rate at low energy. 

This analysis, however, does not apply for the most interesting case of ultra high 
energy photons around 10 19 eV (see below section [8) given that at these energies 
m 2 ^> m 2 if if;' 3 ' > 10~ 17 and ^' 4 ' > 10~ 8 . Hence the above mentioned loop con- 
tributions are at most logarithmic, as the momentum circulating in the fermionic 
loop is much larger than m e . Moreover, in this regime the splitting rate depends 
only on my, the only energy scale present in the problem. One then expects the anal- 
ysis proposed in ||3D to be correct and the splitting time scale to be negligible at 
E y ~ 10 19 eV. 



Lepton pair production 

The process e~ — > e~e~e + is similar to vacuum Cerenkov radiation or helicity de- 
cay, with the final photon replaced by an electron-positron pair. Various combina- 
tions of helicities for the different fermions can be considered individually. If we 
choose the particularly simple case (and the only one we shall consider here) where 
all electrons have the same helicity and the positron has the opposite helicity, then 
the threshold energy will depend on only one LV parameter. In BTl was derived 
the threshold for this reaction, finding that it is a factor ~ 2.5 times higher than that 
for soft vacuum Cerenkov radiation. The rate for the reaction is high as well, hence 
constraints may be imposed using just the value of the threshold. 



7.3.4 LV-modifled threshold reactions: Photon pair-creation 

A process related to photon decay is photon absorption, yy r — > e + e~ . Unlike photon 
decay, this is allowed in Lorentz invariant QED and it plays a crucial role in making 
our universe opaque to gamma rays above tents of TeVs. 

If one of the photons has energy (Oq, the threshold for the reaction occurs in 
a head-on collision with the second photon having the momentum (equivalently 
energy) ku = m 2 /(Oq. For example, if ku = 10 TeV (the typical energy of inverse 
Compton generated photons in some active galactic nuclei) the soft photon threshold 
(Oq is approximately 25 meV, corresponding to a wavelength of 50 microns. 
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In the presence of Lorentz violating dispersion relations the threshold for this 
process is in general altered, and the process can even be forbidden. Moreover, as 
firstly noticed by Kluzniak II 13611 and mentioned before, in some cases there is an 
upper threshold beyond which the process does not occur. Physically, this means that 
at sufficiently high momentum the photon does not carry enough energy to create a 
pair and simultaneously conserve energy and momentum. Note also, that an upper 
threshold can only be found in regions of the parameter space in which the y-decay 
is forbidden, because if a single photon is able to create a pair, then a fortiori two 
interacting photons will do 11311 . 

Let us exploit the above mentioned relation tj e ± = (— )"rj^ between the electron- 
positron coefficients, and assume that on average the initial state is unpolarized. In 
this case, using the energy-momentum conservation, the kinematics equation gov- 
erning pair production is the following ll33l 

^^+l-H>"-'+<->"<^r') 

where E, = <^"' /M"~ 2 and f] = tj'"'/M" _2 are respectively the photon's and elec- 
tron's LV coefficients divided by powers of M, < y < 1 is the fraction of momen- 
tum carried by either the electron or the positron with respect to the momentum k 
of the incoming high-energy photon and (Ot, is the energy of the target photon. In 
general the analysis is rather complicated. In particular it is necessary to sort out 
whether the thresholds are lower or upper ones, and whether they occur with the 
same or different pair momenta 



7.4 Synchrotron radiation 

Synchrotron emission is strongly affected by LV, however for Planck scale LV and 
observed energies, it is a relevant "window" only for dimension four or five LV 
QED. We shall work out here the details of dimension five QED (« = 3) for illustra- 
tive reasons (see e.g. 01371 for the mSME case). 

In both LI and LV cases (33], most of the radiation from an electron of energy E 
is emitted at a critical frequency 

co c = 3 -e B tn (81) 

2 E 

where y(E) = (1 — v 2 (.E)) _1 / 2 , and v(E) is the electron group velocity. 

However, in the LV case, and assuming specifically n = 3, the electron group 
velocity is given by 

= ?!=(\-?L + nV>E- 



"W^n 1 -^" «j- <82) 
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Therefore, v(E) can exceed 1 if 77 > or it can be strictly less than 1 if 77 < 0. This 
introduces a fundamental difference between particles with positive or negative LV 
coefficient 77. 

If 77 is negative the group velocity of the electrons is strictly less than the (low 
energy) speed of light. This implies that, at sufficiently high energy, y{E)- < E/m e , 
for all E. As a consequence, the critical frequency (D~(y,E) is always less than a 
maximal frequency 0) c max 1331 , Then, if synchrotron emission up to some frequency 
©obs is observed, one can deduce that the LV coefficient for the corresponding lep- 
tons cannot be more negative than the value for which a),, max = 0) o bs- Then, if syn- 
chrotron emission up to some maximal frequency <a b s is observed, one can deduce 
that the LV coefficient for the corresponding leptons cannot be more negative than 
the value for which ft)' nax = 0) o b s , leading to the bound (33) 



r? (3) > -— ( 1 ^-— ) ■ (83) 



M /0.34efi\ 3/ 
m e \m e (O obs J 

If 77 is instead positive the leptons can be superluminal. One can show that at 
energies E c > 8 TeV/77 1 ' 3 , y{E) begins to increase fasters than E/m e and reaches 
infinity at a finite energy, which corresponds to the threshold for soft VC emission. 
The critical frequency is thus larger than the LI one and the spectrum shows a char- 
acteristic bump due to the enhanced <a c . 



8 Current Constraints on the QED sector 

Let us now come to a brief review of the present constraints on LV QED and in 
other sectors of the standard model. We shall not spell out the technical details here. 
These can be found in dedicated, recent, reviews such as 01381 . 



8.1 mSME constraints 

It would be cumbersome to summarize here the constraints on the minimal Standard 
Model extension (dimension there and four operators) as many parameters charac- 
terize the full model. A summary can be found in 01281 . One can of course restrict 
the mSME to the rotational invariant subset. In this case the model basically co- 
incides with the Coleman-Glashow one |28l . In this case the constraints are quite 
strong, for example on the QED sector one can easily see that the absence of gamma 
decay up to 50 TeV provides a constraint of order 10~ 16 on the difference between 
the limit speed of photons and electrons ||94"1 . Constraint up to <9(10~ 22 ) can be 
achieved on other mSME parameters for dimension four LV terms via precision 
experiments like Penning traps. 
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8.2 Constraints on QED with 0{E/M) LV 

It is quite remarkable that a single object can nowadays provide the most stringent 
constraints for LV QED with 0(E/M) modified dispersion relations, this object is 
the Crab Nebula (CN). The CN is a source of diffuse radio, optical and X-ray radia- 
tion associated with a Supernova explosion observed in 1054 A.D. Its distance from 
Earth is approximately 1 .9 kpc. A pulsar, presumably a remnant of the explosion, 
is located at the centre of the Nebula. The Nebula emits an extremely broad-band 
spectrum (21 decades in frequency, see 01391 for a comprehensive list of relevant 
observations) that is produced by two major radiation mechanisms. The emission 
from radio to low energy y-rays (E < 1 GeV) is thought to be synchrotron radia- 
tion from relativistic electrons, whereas inverse Compton (IC) scattering by these 
electrons is the favored explanation for the higher energy y-rays. From a theoretical 
point of view, the current understanding of the whole environment is based on the 
model presented in 11401 . which accounts for the general features observed in the 
CN spectrum. 

Recently, a claim of \^ 3 '\ <2x 10~ 7 was made using UV/optical polarisation 
measures from GRBs 01411 . However, the strongest constraint to date comes from 
a local object. In II131II the constraint |^ 3 '| <6x 10~ 10 at 95% Confidence Level 
(CL) was obtained by considering the observed polarization of hard-X rays from the 
CN QH) (see also |fl43l ). 



8.2.1 Synchrotron constraint 

How the synchrotron emission processes at work in the CN would appear in a "LV 
world" has been studied in 1144111391 . There the role of LV in modifying the charac- 
teristics of the Fermi mechanism (which is thought to be responsible for the forma- 
tion of the spectrum of energetic electrons in the CN 01451 ) and the contributions of 
vacuum Cerenkov and helicity decay were investigated for n ~ 3 LV. This procedure 
requires fixing most of the model parameters using radio to soft X-rays observations, 
which are basically unaffected by LV. 

Given the dispersion relations (l48l i and (l49t . clearly only two configurations in 
the LV parameter space are truly different: 77+ • rj_ > and 77+ • r\- < 0, where 77+ 
is assumed to be positive for definiteness. The configuration wherein both r\± are 
negative is the same as the (77+ ■ tj_ > 0, r\ + > 0) case, whereas that whose signs are 
scrambled is equivalent to the case (77+ ■ 77 < 0, T7 + > 0). This is because positron 
coefficients are related to electron coefficients through 77^ = — 77^ 1331 . Examples 
of spectra obtained for the two different cases are shown in Fig.Q] 

A x 2 analysis has been performed to quantify the agreement between models 
and data Q139I . From this analysis, one can conclude that the LV parameters for 
the leptons are both constrained, at 95% CL, to be |?7±| < 10~ 5 , as shown by the 
red vertical lines in Fig. [3] Although the best fit model is not the LI one, a careful 
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Fig. 1 Comparison between observational data, the LI model and a LV one with TJ + • T7_ < 
(left) and r/ + • ?7_ > (right). The values of the LV coefficients, reported in the insets, show the 
salient features of the LV modified spectra. The leptons are injected according to the best fit values 
p = 2.4, E c = 2.5 PeV. The individual contribution of each lepton population is shown. 



statistical analysis (performed with present-day data) shows that it is statistically 
indistinguishable from the LI model at 95% CL 111391 . 



8.2.2 Birefringence constraint 



In the case of the CN a (46 ±10)% degree of linear polarization in the 100 keV — 
1 MeV band has recently been measured by the INTEGRAL mission 111461 11421 . 
This measurement uses all photons within the SPI instrument energy band. How- 
ever the convolution of the instrumental sensitivity to polarization with the detected 
number counts as a function of energy, &(k), is maximized and approximately con- 
stant within a narrower energy band (150 to 300 keV) and falls steeply outside this 
range 111471 . For this reason we shall, conservatively, assume that most polarized 
photons are concentrated in this band. Given afcrab = 1.9 kpc, ki = 300 keV and 
k\ = 150 keV, eq. ( l73l leads to the order-of-magnitude estimate |£ | < 2 x 10~ 9 . A 
more accurate limit follows from (l74t . In the case of the CN there is a robust un- 
derstanding that photons in the range of interest are produced via the synchrotron 
proces, for which the maximum degree of intrinsic linear polarization is about 70% 
(see e.g. 111481 ). Figure |2] illustrates the dependence of 17 on £, (see eq.d74li) for the 
distance of the CN and for 17(0) = 70%. The requirement n(£) > 16% (taking 
account of a 3(7 offset from the best fit value 46%) leads to the constraint (at 99% 
CL) 

|£|<6xl0- 9 . (84) 

It is interesting to notice that X-ray polarization measurements of the CN already 
available in 1978 111491 , set a constraint |£ | < 5.4 x 10~ 6 , only one order of magni- 
tude less stringent than that reported in 111411 . 

Constraint (l84l i can be tightened by exploiting the current astrophysical under- 
standing of the source. The CN is a cloud of relativistic particles and fields powered 
by a rapidly rotating, strongly magnetized neutron star. Both the Hubble Space Tele- 
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Fig. 2 Constraint for the polarization degree. Dependence of 77 on £ for the distance of the CN 
and photons in the 150-300 keV range, for a constant instrumental sensitivity £?(k). 



scope and the Chandra X-ray satellite have imaged the system, revealing a jet and 
torus that clearly identify the neutron star rotation axis [150]. The projection of 
this axis on the sky lies at a position angle of 124.0° ± 0.1° (measured from North 
in anti-clockwise). The neutron star itself emits pulsed radiation at its rotation fre- 
quency of 30 Hz. In the optical band these pulses are superimposed on a fainter 
steady component with a linear polarization degree of 30% and direction precisely 
aligned with that of the rotation axis 11511 . The direction of polarization measured 
by INTEGRAL-SPI in the y-rays is obs = 123° ± 11° (Id error) from the North, 
thus also closely aligned with the jet direction and remarkably consistent with the 
optical observations. 

This compelling (theoretical and observational) evidence allows us to use eq. d76l l. 
Conservatively assuming 0; — o b s = 33° (i.e. 3(7 from 0;, 99% CL), this translates 
into the limit 

|£ (3) |<9xl0- 10 , (85) 

and |^ (3) | < 6 x 10~ 10 for a 2a deviation (95% CL). 

Polarized light from GRBs has also been detected and given their cosmological 
distribution they could be ideal sources for improving the above mentioned con- 
straints from birefringence. Attempts in this sense were done in the past 1193111521 
(but later on the relevant observation [ 1531 appeared controversial) but so far we 
do not have sources for which the polarization is detected and the spectral redshift 
is precisely determined. In 111541 this problem was circumvented by using indirect 
methods (the same used to use GRBs as standard candles) for the estimate of the red- 
shift. This leads to a possibly less robust but striking constraints \£^>\ <2.4x 10 -14 . 

Remakably this constraint was recently further improved by using the INTE- 
GRAL/IBIS observation of the GRB 0412 19A, for which a luminosity distance of 
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Fig. 3 Summary of the con- 
straints on LV QED at order 
0{E/M). The red lines are 
related to the constraints de- 
rived from the detection of 
polarized synchrotron radia- 
tion from the CN as discussed 
in the text. For further ref- 
erence are also shown the 
constraints that can be derived 
from the detection of 80 TeV 
photons from the CN: the 
solid black lines symmetric 
w.r.t. the £ axis are derived 
from the absence of gamma 
decay, the dashed vertical line 
cutting the 7] axis at about 
10~ 3 refers to the limit on 
the vacuum Cerenkov effect 
coming from the inferred 80 
TeV inverse Compton elec- 
trons. The dashed vertical line 
on the negative side of the 
7] axis is showing the first 
synchrotron based constraint 
derived in [144| . 
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85 Mpc (z ~ 0.02) was derived thanks to the determination of the GRB's host galaxy. 
In this case a constraint |^' 3 '| < 1.1 x 10~ 14 was derived £l55|0 



8.2.3 Summary 

Constraints on LV QED 0(E/M) are summarized in Figf3]where also the constraints 
— coming from the observations of up to 80 TeV gamma rays from the CN II 1561 
(which imply no gamma decay for these photons neither vacuum Cherenkov at least 
up to 80 TeV for the electrons producing them via inverse Compton scattering) — 
are plotted for completeness. 



8 The same paper claims also a strong constraint on the parameter §' '. Unfortunately, such a 
claim is based on the erroneous assumption that the EFT order six operators responsible for this 
term imply opposite signs for opposite helicities of the photon. We have instead seen that the CPT 
evenness of the relevant dimension six operators imply a helicity independent dispersion relation 
for the photon (see eq. d54t ). 
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8.3 Constraints on QED with 0(E/M) 2 LV 

Looking back at Table [T] it is easy to realize that casting constraints on dimension 
six LV operators in QED requires accessing energies beyond 10 16 eV. Due to the 
typical radiative processes characterizing electrons and photons it is extremely hard 
to directly access these kind of energies. However, the cosmic rays spectrum does 
extend in this ultra high energy region and it is therefore the main (so far the only) 
channel for probing these kind of extreme U V LV. 

One of the most interesting features related to the physics of Ultra-High-Energy 
Cosmic Rays (UHECRs) is the Greisen-Zatsepin-Kuzmin (GZK) cut off l4lll42l . 
a suppression of the high-energy tail of the UHECR spectrum arising from inter- 
actions with CMB photons, according to py — > A + — > pn°(nn + ). This process has 
a (LI) threshold energy £ t h — 5 x 10 19 (co^/1.3 meV) -1 eV {(Ob is the target pho- 
ton energy). Experimentally, the presence of a suppression of the UHECR flux was 
claimed only recently 12T1I201 . Although the cut off could be also due to the finite 
acceleration power of the UHECR sources, the fact that it occurs at the expected en- 
ergy favors the GZK explanation. The results presented in 015711 seemed to further 
strengthen this hypothesis (but see further discussion below). 

Rather surprisingly, significant limits on B, and tj can be derived by considering 
UHE photons generated as secondary products of the GZK reaction |158lll59 l. This 
can be used to further improve the constraints on dimension 5 LV operators and 
provide a first robust constraint of QED with dimension 6 CPT even LV operators. 

These UHE photons originate because the GZK process leads to the production 
of neutral pions that subsequently decay into photon pairs. These photons are mainly 
absorbed by pair production onto the CMB and radio background. Thus, the frac- 
tion of UHE photons in UHECRs is theoretically predicted to be less than 1% at 
10 19 eV 11601 . Several experiments imposed limits on the presence of photons in 
the UHECR spectrum. In particular, the photon fraction is less than 2.0%, 5.1%, 
31% and 36% (95% C.L) at E = 10, 20, 40, 100 EeV respectively lfT6T1fT62l . 

The point is that pair production is strongly affected by LV. In particular, the 
(lower) threshold energy can be slightly shifted and in general an upper thresh- 
old can be introduced PP . If the upper threshold energy is lower than 10 20 eV, 
then UHE photons are no longer attenuated by the CMB and can reach the Earth, 
constituting a significant fraction of the total UHECR flux and thereby violating 
experimental limits H158lll59l[T63l . 

Moreover, it has been shown [159| that the y-decay process can also imply a 
significant constraint. Indeed, if some UHE photon (Ey ~ 10 19 eV) is detected by 
experiments (and the Pierre Auger Observatory, PAO, will be able to do so in few 
years 11611 ), then y-decay must be forbidden above 10 19 eV. 

In conclusion we show in Fig. 2]the overall picture of the constraints of QED 
dimension 6 LV operators, where the green dotted lines do not correspond to real 
constraints, but to the ones that will be achieved when AUGER will observe, as 
expected, some UHE photon. 

Let us add that the same reasoning can be used to further strength the available 
constraints in dimension 5 LV QED. In this case the absence of relevant UHE photon 
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Fig. 4 LV induced by dimension 6 operators. The LV parameter space is shown. The allowed re- 
gions are shaded grey. Green dotted lines represent values of (77, §) for which the /-decay thresh- 
old £y_ f /«. — 10 19 eV. Solid, blue lines indicate pairs (JJ,£) for which the pair production upper 



threshold k,„ 



10 20 eV. 



flux strengthen by at most two order of magnitude the constraint on the photon 
coefficient while the eventual detection of the expected flux of UHE photons would 
constraint the electron positron coefficients down to \r\^ \ < 10~ 16 (see 11591 11381 
for further details) by limiting the gamma decay process (note however, that in this 
case one cannot exclude that only one photon helicity survives and hence a detailed 
flux reconstruction would be needed). 



9 Other SM sectors constraints 

While QED constraints are up to date the more straightforward from a theoretical 
as well observational point of view, it is possible to cast constraints also on other 
sectors of the SM, most noticeably on the hadronic and neutrino sectors. Let us 
review them very briefly here. 



9.1 Constraints on the hadronic sector 



Being an ultra high energy threshold process, the aforementioned GZK photopion 
production is strongly affected by LV. Several authors have studied the constraints 
implied by the detection of this effect HM HSl ED ES HSU USD- However, a 
detailed LV study of the GZK feature is hard to perform, because of the many astro- 
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physical uncertainties related to the modeling of the propagation and the interactions 
ofUHECRs. 

As a consequence of LV, the mean free path for the GZK reaction is modified. The 
propagated UHECR spectrum can therefore display features, like bumps at specific 
energies, suppression at low energy, recovery at energies above the cutoff, such that 
the observed spectrum cannot be reproduced. Moreover, the emission of Cherenkov 
y-rays and pions in vacuum would lead to sharp suppression of the spectrum above 
the relevant threshold energy. After a detailed statistical analysis of the agreement 
between the observed UHECR spectrum and the theoretically predicted one in the 
presence of LV and assuming pure proton composition, the final constraints implied 
by UHECR physics are (at 99% CL) Q68| 

- 10~ 3 < ri ( p ] < 1(T 6 

-10~ 3 < ri { n 4) < 10- ' (r^ 4) > 0) or < 1(T 6 (tj^ 4) < 0) . (86) 

Of course for dimension five operators much stronger constraints can be achieved 
by a similar analysis (order <9(10~ 14 )). 



9.2 Constraints on the neutrino sector 

LV can affect the speed of neutrinos with respect to light, influence possible thresh- 
old reactions or modified the oscillations between neutrinos flavors. Unfortunately 
we have a wealth of information only about the latter phenomenon which however 
constraints only the differences among LV coefficients of different flavors. In this 
case, the best constraint to date comes from survival of atmospheric muon neutrinos 
observed by the former IceCube detector AMANDA-II in the energy range 100 GeV 
to 10 TeV 111691 , which searched for a generic LV in the neutrino sector 11701 and 
achieved (Ac/c)tj < 2.8 x 10~ 27 at 90% confidence level assuming maximal mixing 
for some of the combinations i, j. Given that IceCube does not distinguish neutrinos 
from antineutrinos, the same constraint applies to the corresponding antiparticles. 
The IceCube detector is expected to improve this constraint to (4c/e)y <9x 10~ 28 
in the next few years 01711 . The lack of sidereal variations in the atmospheric neu- 
trino flux also yields comparable constraints on some combinations of SME param- 
eters 032). 

For what regards the time of flight constraints we have to date only a single event 
to rely on, the supernova SN 1987a. This was a peculiar event which allowed to de- 
tect the almost simultaneous (within a few hours) arrival of electronic antineutrinos 
and photons. Although only few electronic antineutrinos at MeV energies was de- 
tected by the experiments Kamiokall, 1MB and Baksan, it was enough to establish 
a constraint {Ac/c) TOF < 10~ 8 QH or (Ac/c) TOF < 2 x 10~ 9 QH by looking at 
the difference in arrival time between antineutrinos and optical photons over a base- 
line distance of 1.5 x 10 5 ly. Further analyses of the time structure of the neutrino 
signal strengthened this constraint down to ~ 10~ 10 1175111761 . 
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The scarcity of the detected neutrino did not allow the reconstruction of the full 
energy spectrum and of its time evolution in this sense one should probably con- 
sider constraints purely based on the difference in the arrival time with respect to 
photons more conservative and robust. Unfortunately adopting Ac/c < 10~ 8 , the 

SN constraint implies very weak constraints, <jjv < 10 13 and t,\, < 10 34 . 

Threshold reactions also can be used to cast constraints on the neutrinos sector. 
In the literature have been considered several processes most prominently the neu- 
trino Cerekov emission v — ¥ y V, the neutrino splitting v — > V VV and the neutrino 
electron/positron pair production v — » ve~e + . Let us consider for illustration the 
latter process. Neglecting possible LV modification in the electron/positron sector 
(on which we have seen we have already strong constraints) the threshold energy is 
for arbitrary n 

F 2 - 4m * <W 

t th,(n)--? ) V*l) 

0{n) 

withS (B) = § v (£ rt /M)"- 2 . 

The rate of this reaction was firstly computed in fiTI for n = 2 but can be easily 
generated to arbitrary n |48| (see also IU77ID . The generic energy loss time-scale 
then reads (dropping purely numerical factors) 

.. ^co^ /AA 3 *"- 2 ) 

T V-pau-- gAR5 \jj , (88) 

where g is the weak coupling and 0„. is Weinberg's angle. 

The observation of upward-going atmospheric neutrinos up to 400 TeV by the 
experiment IceCube implies that the free path of these particles is at least longer 

(3) 
than the Earth radius implies a constraint Tfy < 30. No effective constraint can be 

optioned for n = 4 LV, however in this case neutrino splitting (which has the further 

advantage to be purely dependent on LV on the neutrinos sector) could be used on 

the "cosmogenic" neutrino flux. This is supposedly created via the decay of charged 

pions produced by the aforementioned GZK effect. The neutrino splitting should 

modify the spectrum of the ultra high energy neutrinos by suppressing the flux at 

the highest energies and enhancing it at the lowest ones. In [178] it was shown 

that future experiments like ARIANNA 01791 will achieve the required sensitivity 

to cast a constraint of order r\\, ' < 10~ 4 . Note however, that the rate for neutrino 

splitting computed in 11781 was recently recognized to be underestimated by a factor 

0(E/M) 2 [ 180|. Hence the future constraints here mentioned should be recomputed 

and one should be able to strengthened them by few orders of magnitude. 



10 Summary and Perspectives 

We can summarize the current status of the constraints for the LV SME in the fol- 
lowing table. 
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Table 2 Summary of typical strengths of the available constrains on the SME at different orders. 

Order I photon \e~/e + iHadrons I Neutrinos" 



n=2 
n=3 
n=4 



N.A. 

o(irr 14 ) (GRB) 

0(l(r 8 ) (CR) 




0(1(T 8 ) 
O(30) 
0(1(T 4 )* (CR) 



GRB=gamma rays burst, CR=cosmic rays 

" From neutrino oscillations we have constraints on the difference of LV coefficients of different 
flavors up to 0(1O~ 28 ) on dim 4, 0(1O~ 8 ) and expected up to 0(1O~ 14 ) on dim 5 (ICE3), expected 
up to 0(1O~ 4 ) on dim 6 op. * Expected constraint from future experiments. 



A special caveat it is due in the case of n = 4 constraints. As we have seen, they 
mostly rely (in the QED and Hadronic sector) on the actual detection of the GZK 
feature of the UHECR spectrum. More specifically, UHECR constraints have relied 
so far on the hypothesis, not in contrast with any previous experimental evidence, 
that protons constituted the majority of UHECRs above 10 19 eV. Recent PAO 111811 
and Yakutsk II 1821 observations, however, showed hints of an increase of the aver- 
age mass composition with rising energies up to E rj 10 196 eV, although still with 
large uncertainties mainly due to the proton-air cross-section at ultra high energies. 
Hence, experimental data suggests that heavy nuclei can possibly account for a sub- 
stantial fraction of UHECR arriving on Earth. 

Furthermore the evidence for correlations between UEHCR events and their po- 
tential extragalactic sources 01571 — such as active galactic nuclei (mainly Blasars) 
— has not improved with increasing statistics. This might be interpreted as a further 
hint that a relevant part of the flux at very high enrages should be accounted for 
by heavy ions (mainly iron) which are much more deviated by the extra and inter 
galactic magnetic fields due to their larger charge with respect to protons (an effect 
partially compensated by their shorter mean free path at very high energies). 

If consequently one conservatively decides to momentarily suspend his/her judg- 
ment about the evidence for a GZK feature, then he/she would lose the constraints 
at n = 4 on the QED sectotjj as well as very much weaken the constraints on the 
hadronic one. 

Assuming that current hints for a heavy composition at energies E ~ 10 196 eV 
111811 may be confirmed in the future, that some UHECR is observed up to E ~ 
10 20 eV 111841 . and that the energy and momentum of the nucleus are the sum of 
energies and momenta of its constituents (so that the parameter in the modified 
dispersion relation of the nuclei is the same of the elementary nucleons, specifically 



9 This is a somewhat harsh statement given that it was shown in 11831 that a substantial (albeit 
reduced) high energy gamma ray flux is still expected also in the case of mixed composition, so 
that in principle the previously discussed line of reasoning based on the absence of upper threshold 
for UHE gamma rays might still work. 
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rj p ) one could place a first constraint on the absence of spontaneous decay for nuclei 
which could not spontaneously decay without LV| 10 I 

It will place a limit on r\ p < 0, because in this case the energy of the emitted 
nucleon is lowered with respect to the LI case until it "compensates" the binding 
energy of the nucleons in the initial nucleus in the energy-momentum conservation. 
An upper limit for rj p > can instead be obtained from the absence of vacuum 
Cherenkov emission. If UHECR are mainly iron at the highest energies the con- 
straint is given by r\ p < 2 x 10 2 for nuclei observed at 10 196 eV (and r\ p < 4 for 
10 20 eV), while for helium it is 7] p < 4 X 1(T 3 lfl85l 

So, in conclusion, we can see that the while much has been done still plenty 
is to be explored. In particular, all of our constraints on 0(E/M) 2 LV EFT (the 
most interesting order from a theoretical point of view) are based on the GZK effect 
(more or less directly) whose detection is still uncertain. It would be nice to be able 
to cast comparable constrains using more reliable observations, but at the moment 
it is unclear what reaction could play this role. Similarly, new ideas like the one of 
gravitational confinement ||89l presented in section 14.1.21 seems to call for much 
deeper investigation of LV phenomenology in the purely gravitational sector. 

We have gone along way into exploring the possible phenomenology of Lorentz 
breaking physics and pushed well beyond expectations the tests of this fundamental 
symmetry of Nature, however still much seems to await along the path. 

Acknowledgements I wish to that Luca Maccione and David Mattingly for useful insights, dis- 
cussions and feedback on the manuscript preparation. 
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